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We consider geometrical characteristics of monopole clusters of the lattice SU{2) 
gluodynamics. We argue that the polymer approach to the field theory is an ade¬ 
quate means to describe the monopole clusters. Both finite-size and the infinite, or 
percolating clusters are considered. We find out that the percolation theory allows 
to reproduce the observed distribution of the finite-size clusters in their length and 
radius. Geometrical characteristics of the percolating cluster reflect, in turn, the 
basic properties of the ground state of a system with a gap. 
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I. INTRODUCTION 

Explanation of the conhnement in terms of the monopole condensation was proposed as 
early as aronnd the year 1974 [^. Moreover, the idea is strongly snpported by the lattice 
data 1^. Nevertheless, nnderstanding the lattice data in terms of the continunm theory 
still represents a challenge. Indeed, generically one nsnally thinks in terms of a Higgs-type 
model: 

= j d‘x{lF'^, + + V(\4,\^)) (1) 

where 0 is a scalar held with a non-zero magnetic charge, is the held strength tensor 
constrncted on the dnal-gluon held B^, is the covariant derivative with respect to the 
dnal gluon. Finally, U(|0p) is the potential energy ensuring that (0) 7 ^ 0 in the vacuum. 
There exist detailed hts to the numerical lattice data within such a framework |^. However, 
relation of the “ehective” helds (f, to the fundamental QCD helds remains unclear. 

One of the main problems is the lack of understanding of the nature of the monopole held 
conhgurations. In particular, there is seemingly no reason to think about the monopoles as 
quasi-classical objects, for review see @]. Under the circumstances, it is natural to exploit 
the lattice measurements to accumulate hrst data on, so to say, phenomenology of the 
monopoles. In fact, dehnition of the monopoles in the non-Abelian case is not unique and 
we will concentrate on the monopoles dehned within the Maximal Abelian Projection (MAP) 
in SU(2) lattice gauge model. 

A gauge-variant dehnition of the monopoles would not be of much interest by itself. 
The basic idea is that there are gauge-invariant objects behind which are detected through 
the projection. This, gauge-invariant facet is manifested in gauge invariant properties of 
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the MAP monopoles. In particular, (see and references therein) the three dimensional 
monopole density Pmon does not depend on the lattice spacing and is given in the physical 
units: 

Pmon = 0.65(2) ~ ’ (2) 

where (Jsu{ 2 ) is the string tension. 

An important remark is in order here. While discussing the monopole density one should 
distinguish between hnite-size clusters and the percolating cluster P, |^. There is a spectrum 
of the hnite-size clusters, as a function of their length, while the percolating cluster is in a 
single copy. In other words, the percolating cluster hlls in the whole of the lattice and its 
length is proportional to the volume of the lattice V 4 : 

Lperc Pmon ' ^4 (3) 

The observation (^ refers only to the percolating cluster. 

Also, upon identihcation of the monopoles in the Abelian projection, one can measure 
the non-Ahelian action associated with these monopoles. The results [Q turned in fact 
astonishing and can be explained only as a hne-tuning |^. To explain the meaning of the 
hne tuning, let us remind the reader that the probability of hnding any held conhguration 
is a product of action- and entropy- factors. Then, it turns out that for the monopoles at 
the presently available lattices both factors diverge in the ultraviolet but cancel each other 
to a great extent: 


Pmon = exp(-S'mon) X (entropy) 


exp(—Cl ■ L/a) ■ exp(-(-C2 ■ L/a), 

|ci - C2I 


A, 


QCD 


( 4 ) 

( 5 ) 


where Pmon is the probability to hnd the monopole trajectory of the length L, a is the lattice 
spacing and ci ,2 could in principle depend on the a as well but are treated as constants to 
describe the data on the available lattices (with smallest lattice spacing a ~ 0.06 fm from 
Ret. 0). 

The observation (^) implies also that the monopoles look point-like on the presently 
available lattices. Indeed, let us normalize the phenomenological monopole action for the 
lattice monopoles to the case of point-like monopoles. For the Dirac monopole the action 
would be: 

CXD 

S = L ■ — [ = const. , ( 6 ) 

Stt J ae^ 


where e is the electric charge and B is the magnetic held of the monopole, ~ 1/e^r"^. 

Thus, at presently available lattices the MAP monopoles are structure-less, or point-like 
objects. This observation immediately triggers further questions since usually one considers 
generic non-perturbative huctuations as large-scale. We would not discuss such questions 
now but simply stick to the phenomenological observation that at the presently available 
lattices monopoles can be treated as geometrical objects, with no size. One can advocate 
then P] the use of the polymer approach to the held theory, see, e. g.. Ref. m- 

Our basic idea here is to pursue a geometrical language to describe monopoles, without 
any explicit use of Lagrangians like (|I[) . 

In Sect. H we consider the relation between the percolation theory and the monopole 
physics. In Sect. |IT1| we use the percolation theory to evaluate the spectrum of the hnite- 
size monopole clusters. In Sect. E the same hnite-size clusters are described in terms of 
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the simplest vacuum loop in the polymer representation. In Sect.0 we comment on the 
properties of the percolating cluster as representing the ground state of the system. In 
Sect.|^ we discuss measurements which could clarify further the structure of the effective 
theory describing the monopoles. In Sect. El we present conclusions. 

It is worth mentioning that in a few cases we reproduce in some detail argumentation 
well known in other helds, in particular in the condense matter. And in this sense our 
presentation is not original in these cases. However, we believe that the justihcation of the 
use and application of the condense-matter techniques to non-perturbative fluctuations of 
the Yang-Mills helds is new. 


II. MONOPOLE CLUSTERS AND PERCOLATION 


In this section we will outline interconnections between properties of monopole clusters, 
percolation theory [|^ and held theory. The aim is to motivate our basic assumptions (see 
subsection EB and to apply later the percolation theory to the monopole clusters. 


A. Monopole condensation in compact 17(1) 


A two-line theory of the monopole condensation was presented hrst in Ref. for the 
case of the compact (or lattice) U{1) theory. We will review the derivation here and dwell 
on its connection with percolation. 

The role of the compactness of the 17(1) lattice theory is to ensure that the Dirac string 
does not cost any energy (for a review see, e.g., Ref. [Q). That is why in Eq. (H) we take 
into account only the energy of the radial held. Nevertheless, the monopoles are inhnitely 
heavy in the limit a —> 0 and, at hrst sight, this precludes condensation since the probability 
to hnd a monopole trajectory of the length L is suppressed as 


exp(—S') = exp ( — 


a 


(7) 


Note that the constant c depends on details of the lattice regularization but can be found 
explicitly in any particular case. 

What makes the condensation feasible, is an exponentially large enhancement factor due 
to the entropy. Namely, a trajectory of length L for a point-like monopole can be realized 
on a cubic lattice in Nl = various ways. To evaluate the Yl let us notice that the 
monopole occupies the center of a cube and at each step the trajectory can be continued to 
an adjacent cube. In four dimensions there are 8 adjacent cubes. However, one of them is 
to be excluded from the counting since the monopole world-line is non-backtracking^. Thus 
the entropy factor is: 


Nl 


exp 



( 8 ) 


^ If a piece of the trajectory is covered in the both directions it is not observed on the lattice. Physically, 
this cancellation corresponds to the cancellation between monopole and anti-monopole. 
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and it cancels the snppression dne to the action (|^ if the conpling satishes the condition 
of criticality: 

eLi = c/ln7 ^ 1 (9) 


where we qnote the nnmerical valne of for the Wilson action and cnbic lattice. At et 


■'crit ^ OXJ.V... V V IIOV/IX CXV.^UlV/iX CXiXVJ. XCXU X XU 

any trajectory length L is allowed and the monopoles condense. This simple framework is 
verihed within abont one percent accnracy as far as the prediction of is concerned 


13, 


14 


One can say that the conpling of the compact 17(1) is to be hne-tnned to trigger the 
phase transition. 


B. Relation to percolation 


In fact the derivation of the preceding snbsection can be viewed as an application of the 
percolation theory |]^. Moreover, one thinks in terms of the simplest percolation possible, 
that is, nncorrelated percolation. 

Indeed, monopoles are observed as trajectories on the links of the dnal lattice |^. Pos- 
tnlate that the probability to occnpy a link is given by: 


p = exp(—c/e^) 


( 10 ) 


compare Eq. at L = a. The probability that (nncorrelated) links form a trajectory of 
length L is then given by Eq. ( 0 )- 

Formation of an inhnitely long, or percolating clnster at a critical valne p = pcrit is a 
common featnre of percolating systems. In onr case. 


Vcrit = 1/7 , (11) 

see Eq. (||). 

It might worth mentioning that in text-books one wonld rather End Pcrit = 1/8 since 
the non-backtracking featnre of the monopole trajectory is specihc for charged particles. 
From the theoretical point of view the non-backtracking is a manifestation of the monopole 
charge conservation. Another conseqnence of this conservation law is the closeness of the 
trajectories which has not been taken into acconnt so far. As we shall see in Sect. |V| the 
closeness of the trajectories brings in only a pre-exponential factor and can be ignored at 
the moment for this reason. 


C. Relation to field theory 


The derivation in Section [11 A| implies that the monopole condensation occnrs when the 
monopole action is nltraviolet divergent. On the other hand, the onset of the condensation in 
the standard held theoretical langnage corresponds to zero mass of the magnetically charged 
held 0. It is important to realize that this apparent mismatch between the two langnages 
is not specihc for the monopoles at all. Actnally, there is a general kinematic relation (see, 
e.g.. Ref. [jl^ 1^) between the physical mass of a scalar held and the mass M defined 

in terms of the (Enclidean) action S, M = S/L: 





8 . 


^Irop-a ^ C‘i-{M{a) 


a 


( 12 ) 
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It is rni that enters the propagator of a scalar particle DijP'), 


DsciP^) 


1 

’ 


(13) 


not M‘^{a). The factor In 7 in Eq. (^) is specihc for the “monopole walk” (see the preceding 
subsection). The Eq. ([T^) establishes a connection between uncorrelated percolation and 
propagator of a free particle (in Euclidean space), 


~ MPcrit/p), Pcrit>P, 


where we neglect, as usual, higher corrections in a. 

It is also worth mentioning that the overall normalization of the propagator of a free 
particle in the polymer approach reads in the leading (in the lattice units a) order as follows: 


exp|-M(a)Lp„t,j| = Cl^a? Dsc{xi,Xf ), (14) 

paths 


where Dscaiar{x,x) is the standard held-theoretical propagator of a scalar particle with 
mass (|l^ in the Euclidean space-time, Lpath is the length of a particular path connecting 
the initial and hnal points Xjj . Note that the action factor for a given path, exp(—M -Lpai/i), 
is similar to the action of a polymer of length L, with the chemical potential replaced by 
the mass M(a) 

Eq. ([T^) demonstrates that the propagating mass can be kept hnite in the limit a —> 0 
only at price of hne tuning. Moreover, Eq. (0) looks exactly the same as the standard 
radiative correction to the mass of a scalar particle. In more detail, the hrst term in the 
right-hand side is actually the magnetic held energy. The In 7 term, on the other hand, plays 
the role of a counter-term and this counter-term is calculable in our case! 


D. Percolation and the quadratic divergence in field theory 


Now we come to explain an important difference between phase transitions in case of 
uncorrelated percolation and in case of the monopoles. This difference might come as a kind 
of surprise since the transition point can be determined within the percolation theory with 
a good accuracy (see the discussion of the U{1) case above). 

The difference concerns properties of a percolating system and of the monopole vacuum in 
case of p > Per and in the conhning phase, respectively. If p > Per the product of probability 
and entropy factors reduces to exp(-(-|£| ■ L), where e can be arbitrarily small. Then, it is 
clear that the longest possible trajectory wins and the cluster Ells in a hnite fraction of the 
whole lattice. In the percolation-theory language the statement is that the fractal dimension 
at p > Perit coincides with the space dimension D: 


i,Dfr')perc ^ i P ^ Perit j 


(16) 


see, e.ff., Ref. 0 and references therein. Note that in the language of held theory p > pern 
corresponds to a tachyonic mass, < 0, see Eq. m- 

In case of monopoles, Eq ([I^ ) would correspond to 


i^Pmon) 


,-3 


mon Jperc 


( 16 ) 
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Note that in the held theoretical language Eq. (|T^) would be interpreted as a power-like 
divergence well known in perturbation theory of charged scalar particles. Indeed, already 
on dimensional considerations one would conclude that: 

(|*/*| )perc Oj^Pmon)perc ~ ® 5 

(for a more careful derivation see Ref. P). 

The behaviour ( p!6[) is in a sharp contrast with experimental data, see and this is 
what we mean by the difference in the behaviour of a percolating system at p > Pcr and 
of the monopoles in the conhning phase. In other words, the hne tuning exhibited by the 
monopoles in the conhning phase is a specihc feature of the vacuum state of the lattice gauge 
theory. 


E. Fine tuning and Coulomb-like interaction 

Let us emphasize again that in the present note we treat the hne tuning of the monopoles 
as a pure observation and look for its implications. We are not in position at all to track 
its origin back to the non-Abelian Lagrangian which after all determines the results of all 
the lattice measurements. One of very few theoretical points which we can nevertheless add 
is that a long-range. Coulomb-like force is a necessary ingredient of a hne-tuned theory. Of 
course, in case of color by long distances we mean r S> a but r < 

To substantiate the point, let us consider the conhning phase of the compact U{1). 
The relative simplicity of this case is that the dynamics is explicitly known and, moreover, 
determined by the value of the electric charge e which can be tuned arbitrarily “by hand” 
(while in the non-Abelian case the coupling runs and is not under control in this sense). 

Thus, one can introduce two scales in the conhning phase of the compact U{1) in the 
following way: 

Ruv = 0 , , Rjr = — , e = (e^ — > 0 ; e -Cl. (18) 

At hrst sight, we are in the same situation as in case of uncorrelated percolation since we 
introduced a tachyonic mass, 

e ^ e ^ , /i~-, 

a 

and nothing can prevent development of a quadratic divergence, see discussion above. 

On the other hand it is known from the lattice measurements (see, e.g., Ref. and 
references therein) that the IR scale does coexist with the UV one. And, indeed, there is 
the loophole in our logic that introducing the tachyonic mass would immediately result in 
a quadratic divergence. This would be true if the interaction were local. However, imagine 
that there are other monopoles at a distance of order Rjr. Then these monopoles can modify 
the Coulomb field of the monopole considered by order unit at r ~ Rir. This, in turn, can 
be interpreted as a change in the mass (due to interactions): 

\6M(a)\ ~ -. (19) 

a 

This change can compensate for the “short-distance” tachyonic mass so that the particles 
are kept away from each other at a distance of order Rir. Which also means that there is 
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no power-like UV divergences (see discussion above). This double-face interpretation of (|T^ 
as IR and UV effects is a unique feature of the Coulomb-like interaction. 

Note that we are not really giving a proof that the introduction of Rjr is self-consistent. 
We just argue that without long-range forces the hne tuning would be not possible at all. In 
case of non-Abelian theories there are many more open questions since the monopole action 
is not bounded from below (see, e.g., Ref. |^) and the fact that the action is tuned to the 
entropy is even more difficult to explain theoretically. 


F. Fine tuning and models 


It is worth to emphasize that the hne tuning which is observed for the lattice monopoles 
is of the same generic type considered so mystifying in case of the Standard Model. There 
is, however, a peculiarity in our case. Indeed, if we compare the “natural” estimate for Pmon 
( 0 ) with the data (^ we see that there are three powers of a ^ that are “tuned away”, not 
two powers as in the standard held-theoretical language. And, indeed, the data (^j imply 
in the “naive” a —>■ 0 limit [Q: 

lim{|(/)|^) ~ Pmon -a ^0, 

a—>0 

limm^ ~ lim-(M(a) - In 7/a) ~ ^QCD ^ 
a^o a^o a a 

a—>^0 ^ 


where nigiuon is the (dual) gluon mass which arises in theories of the type (|l|). 

Eqs. (^^ are of course not what one would expect for the standard theory with 
spontaneous symmetry breaking. It is worth emphasizing, however, that the limit a —>■ 0 in 
Eqs. (^^ should be understood with some reservations. What we mean actually in Eqs. (^0[) 
by a —> 0 limit is “the lattice spacing a as small as possible within availability on present 
lattices”. The behaviour (^0[) can actually change in the academic limit a —0 (i.e. at 
lattice spacings much smaller than those presently available). Subsequent considerations in 
the present paper would actually suggest such a possibility (see Sect. 0 ). 

Although Eqs. (^) may not survive at smaller a, they do answer questions why new 
point-like fluctuations (implied by the hne-tuning) do not disturb the standard /3-function 
of the SU{2) gluodynamics at existing lattices. Indeed, according to (^OD all the scalar 
degrees of freedom are actually removed from the physical spectrum if a —> 0. 

The properties (pOD imply also that the potential energy U(|(/)p) scales in the limit a —>■ 
and our effective theory can well be a useful approximation to study the vacuum properties. 
In particular, the monopole conhning mechanism survives in the limit |P[ a ^ 0. 


G. Formulating the main hypothesis 

After all these preliminary discussions we are set to formulate our main hypothesis. 
Namely, we will assume that we can consider the point 

P = Pcrit (21) 

as adequately describing the physics in the confining phase of the non-Abelian gauge theory. 
The justihcation for this hypothesis is the hne tuning observed on the lattice and discussed 
in length above. 






There are important reservations to be made. Namely, the fine tuning does imply that 
the physics is so to say “frozen” at p = Pcr as far as the UV scale is concerned. However, as 
far as the dependence on the scale Aqqd is concerned it could be different than at the point 
m- Moreover, we can give a more quantitative meaning to the scale ^''Aqcd" in fhe case 
considered. The point is that the percolating cluster has self-crossings and the length of the 
trajectory between these self-crossings can be considered as the monopole free-path length. 
Direct measurements show [|^ that this length (measured along the trajectory) scales: 

Lfree ~ 1.6/m . (22) 

Thus, we can apply our hypothesis as long as 

Lfree '> L a. (23) 

In the next section we exploit the percolation theory to describe the structure of the 
finite-size monopole clusters satisfying (|23|) . 


III. FINITE-SIZE CLUSTERS AND PERCOLATION 


A. Data and percolation picture 


Detailed data on the structure of the monopole clusters were obtained in 0]. As was 
mentioned above, there is a single percolating cluster, whose length grows with the lattice 
volume, and finite-size clusters. In this section we will concentrate on the finite-size clusters 
satisfying the condition (^). These clusters are characterized, first of all, by their length. 
It was found that the length spectrum is described by a power law: 

N(L) = (24) 

where 

T ^ 3 (25) 

for all lattice spacings and sizes tested and the coefficient C 4 depends only weakly on /?. For 
our purposes we can neglect this dependence. 

Another important characteristics of the clusters is their radius Rl as function of the 
length L. By the radius one understands the average distance between two cluster links: 


Ri = 


a 


Lja 




a 


If a 


i=l 


Yl N- 


(26) 


where Xjj are coordinates of the links and x = ( 1 /T) The measurements indicate: 

Rl ~ const 1 -|-const 2 V^. (27) 

Thus, our problem is to clarify whether the data (p5|) and (p7| ) can be understood within 
the percolation theory (and our main hypothesis, see subsection |11 h] ). It is encouraging to 
observe that even without any dynamical input we can conclude that the lattice data on 
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the finite-size monopole clusters reveal a picture typical for percolating systems. Indeed, a 
generic form of the spectrum, for p < p^rit is 

N(L) ~ ( 28 ) 

where r is the so called Fischer exponent and p vanishes at the critical point, p = Pcrit- The 
latter is easy to understand since at the point of the phase transition the correlation length 
is infinite and there is no dimensional parameter left. Moreover, a power law like (|?^) is 
common to the percolating systems, 

Rl ~ ( 29 ) 

and Dfr is called the fractal dimension. For the monopoles we apparently have from the 
lattice measurements: 

{Djr)mon ~ 2 . ( 30 ) 

Thus, the data on the finite-size monopole clusters exhibit a typical percolation picture 
and the next question is whether it is possible to evaluate the exponents (^3]) and (|3(]|). 


B. Hyperscaling relation 


Similarity between percolation and properties of dynamical systems undergoing the phase 
transition is well known []T|. However, direct evaluation of the critical exponents requires, 
generally speaking, a particular dynamical input on the system considered. Mostly, knowl¬ 
edge on the excitation spectrum is required. 

However, there is a general relation between the critical exponents which we are going to 
apply first. Note that our considerations here are not specific for the monopoles and similar, 

[13 


e.g., to the analysis of Ref. 


Still, for the sake of orientation let us mention a few points 


important for the application of the percolation theory to the monopoles clusters. 

First, one derives all relations for p < p^r, without considering p > pcr- As we discussed in 
length in Sect. |I|, uncorrelated percolation can be indeed relevant to the monopole physics 
only at p < pcr- Second, all the relations of the percolation theory follow more or less directly 
from the assumption ( p8[) and it is worth emphasizing that the spectrum (^) is very natural 
physics-wise. Indeed, the presence of the exponential factor at p — pcr < 0 is obvious from 
considerations of subsection [11 H|. Moreover at the point p = pcr there is no scale left and 


the power like dependence of the spectrum on L is the only dependence allowed. In this way 
one comes to include the factor L~'^ as well. 

Starting from the spectrum (P^) it is quite straightforward, see, e.g., Ref. |^, to derive 
the following well-known relation: 


D 


T = 


D 


+ 1 


fr 


(31) 


where d is the dimension of space-time and Dfj. is in fact Dfr{p = Pcr)- 
In our case, 

T ^ 3 , D = A, Dfr ^ 2 . 

Thus, Eq. (1^) is satisfied within the error bars of the lattice measurements and this 


(32) 


observation is one of our main results. In view of its importance, we will later rederive 
(P^) and some generalizations of it in the language of field theory. 
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C. Fractal dimension 


As is mentioned above, the relation between the length and radius of the cluster is deter¬ 
mined by the fractal dimension, see Eq. (^Tl). The fractal dimension, in turn, is determined 
by the kind of the walk. 

In fact we are dealing with the monopole walk which, to the best of our knowledge, has 
not been studied in detail in the literature. However, the characteristics of the monopole 
walk are so to say flanked by the characteristics of the well known random and self-avoiding 
walks. Indeed, the monopole walk chooses freely one of 7 directions available at each step. 
This is in common with the random walk. On the other hand, choosing the eighth direction 
would result in an immediate self-crossing and this is forbidden. The latter feature is in 
common with the self-avoiding walk. However, in contradistinction from the self-avoiding 
walk, self-crossings are allowed for the monopole trajectories at later stages. 

The observation central for this section is that in 77 = 4 the fractal dimension = 2 
both for the random and self-avoiding walks. Therefore we can predict Df^ = 2 for the 
monopole walk as well. 

In more detail, for the random walk one has in any number of space dimensions 


fr)random 2 . 

For the self-avoiding walk one has (the so called Flory’s fractal dimension): 


(33) 


(34) 


which is valid for D < 4 and gives exactly the same D= 2 at D = 4. 

Thus, the experimental value {Dfr)man = 2 appears to be well understood theoretically. 


IV. SEEING FREE MONOPOLES AT SHORT DISTANCES 

From our discussion of the relation between the uncorrelated percolation and free-held 
theory (see Sect. Eig) it is clear that the results of the preceding section can be derived 
in terms of Feynman graphs as well Moreover, since we are considering relatively short 
trajectories, see (p3D , the effect of the mass can be neglected. In this section we will reproduce 
the results above by considering the simplest Feynman graph relevant, that is a single 
vacuum loop. This approach allows us to derive also some generalizations, like length-of- 
the-trajectory spectrum at hnite temperature. 


A. Polymer representation for a massless particle propagator 


In this subsection we will evaluate the Fischer exponent r (see (^^) starting from the 
observation that the huctuations of the held 0 are massless on the scale a. To this end we 
reproduce hrst the basic points of the polymer approach to the held theory of a free scalar 
particle, see, e.g., Refs. |]^, |^, ^ . 


^ This relation, in connection with the monopole clusters, is in fact mentioned in Ref. |Q . 

^ We are grateful to D. Diakonov for a discussion on this point and providing us with a reference to pH|. 
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The partition function for a closed polymer is: 


Z = 


/ CXD 

7V=1 


-M-N 


N 


Zn{x,x) , 


(35) 


where M is the chemical potential and Z]s[{xo,Xf) is the partition function of a polymer 
broken into N segments: 


N-l „ AT p 

ZN{xo,Xf) = [JJ / JJ 

i=l i=l - 


(5(| xi - Xi-i \ - a) 
27r^a^ 


(36) 


This partition function (^5]) contains a summation over all atoms of the polymer weighted 
by the Boltzmann factors. The ^-functions in (^) ensure that each bond in the polymer 
has length a. The starting point of the polymer is Xq and the ending point is Xf = xp^. 

Note that there is a pre-exponential factor 1/N in Eq. (|35D. This is due to consider¬ 
ing closed trajectories. Indeed, the factor is introduced to compensate for the N-multiple 
counting of the same closed trajectory in the partition function (|35D since any atom on this 
trajectory can be considered as the initial and hnal point. As we shall see later, our hnal 
result crucially depends on the pre-exponential factors. 

The crucial step to relate (^) to a free particle path integral is the so called coarse- 
graining. Namely, the A^-sized polymer is divided into m units by n atoms {N = mn), and 
the limit is considered when both m and n are large while a and ^/na are small. We get. 


{v+l)n—l 


n 


27r^a^ 


6{\xi - Xi+i\ - a) ^ f— 

V Trna^ / 


exp 


na‘ 

,th 




X,. 


(37) 


function becomes [22 


where the index i, i = un - ■ ■ + — labels the atoms in unit. The polymer partition 

{xy - x^_i)2- 


Z]\f{xo,Xf) = const- 


m—1 

Ifd' 

U=1 


X 


, 9 . 2m , rn 

(^) 


U=1 


na^ 


■ exp 


n ■ a n 

i/=i 


(38) 


The Xj’s have been re-labeled so that x^y is the average value of x in at the coarser cell. 
Note also that at this stage there appears the chemical potential /i related to the original 
parameter M through 

/i = M —In 7/a, (39) 

as is discussed above. 

Using the variables: 


1 2 7^2 i\T 2 8p 

s = - na u, 1 = -a N, = — 


(40) 


one can rewrite the partition function (^^ as 


dl 


Z = const ■ / — 


I Dx exp<( — 

0 x(0)=x(l)=x 


^xj(s) + m. 


2 

prop 


ds 


(41) 
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Note that the mass renormalization in Eq. (^OD is consistent with Eqs. ([T^j39|) . 

After these preliminary steps we can readily derive the distribution in the length of the 
trajectories in the massless case. To this end, let us rescale and s in such a way, that 
there is no I dependence left in the action if niprop = 0: 


L = l/a, s = s/l, Xfj, 


x^/vi . 


( 42 ) 


Then, indeed. 


where 


Z 


oo 


const 


dL 

T 




I = J Dx exp 

:r(0)=x(/)=i 


-\J ^liS) ds 
0 


( 43 ) 


( 44 ) 


At hrst sight, Eq. (^B] ) implies that we have a.dL/L spectrum since there is no L dependence 
left otherwise. However, the actual spectrum refers to the number of loops in a unit volume 
(see and the discussion of it). Thus, the dL/L spectrum refers to the volume in the x 
units. Since the x^ ~ we have in the physical-volume units: 


N(L)dL = CO Vo , ( 45 ) 

where Vn is the volume in D-dimensional space and cu is a constant. Although we are 
interested only in the D = A case we kept D as variable to emphasize that we rederive in 
fact the hyperscaling relation m- Note that Djr = 2 is implicit in our derivation and is 
encoded in fact in the transformation (|^ . 


B. Coulomb-like interaction 


So far we considered approximation of free particles which corresponds to the dominance 
of the monopole self-energy. We expect, however, that the monopoles interact also Coulomb¬ 
like. Other, effective interactions are not ruled out either. The Coulomb-like interaction can 
readily be included into the action in the polymer representation. The corresponding extra 
piece in the action is given by: 

Scoulomb „ f f dS\dS2X\^p,Dpy{X\ 3)2)h2,i/ (^6) 

^ Jo Jo 


where qm is the magnetic charge, Dp^{z) is the massless photon propagator, and the dot, 
as usual, means differentiation with respect to the proper time. 

The action (^) is manifestly invariant under the rescaling (^): 


5, 


Coulomb 


9m 



dsidhxi,pD 


X2)X2,V 


0 Jo 


(47) 
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and, therefore, the 1 /behaviour of the spectrum (^51) should be still true for^ D = A upon 
inclusion of this interaction^. 


C. Finite temperatures, finite volume lattices 


Formalism of Section [IV A| can readily be generalized to the case of a non-zero tempera¬ 
ture. Indeed, we should rewrite now in terms of a length-of-the-trajectory distribution the 
propagator of a massless particle at hnite temperature. 

As usual, hnite temperature in Euclidean space corresponds to a compactihed fourth 
direction. Thus, we simply write the L-distribution corresponding to an ensemble of non¬ 
interacting particles with masses = (27rn)^T^ in d=3 space-time: 


N{L) = C3 14 -pyj J]] exp ( - {2TTnTfL ■ a), 


(48) 




where V3 = V4T is the three dimensional volume of the time-slice. Note that in the limit 
T —>■ 0 the sum over the exponentials is proportional to 1/ (Ta/L ■ a) and we come back to 
(^) upon proper normalization of the constant C 3 , 


C3 = 2 C4 v 7 ra , 


where C4 enters Eq. (|24D. However, in the opposite limit, T —> 00 we effectively get three 
dimensional theory with the volume V3 and the loop distribution N{L) oc 

A remark on the range of validity of (]48|) is now in order. The point is that in (|48|) we 
disregard chemical potential p 7^ 0. Whether this is a valid approximation depends on the 
numerical values of the parameters involved, /i, T, a. In particular, if we tend a ^ 0 and 


keep p in physical units the approximation (|48|) is not valid. 

Proceeding in a similar way we can derive also the spectrum on a hnite-size lattice of the 
volume V4 = n^=i ^1 
periodic boundary condition we get: 


where is the length of the space in direction. Assuming the 


/ 4=1 


exp { - ■ a} , 


(49) 


where 


Co = C4a^ . 


V. INFRARED CLUSTER AND PROPERTIES OF THE GROUND STATE 


In this section we discuss geometrical properties of the percolating cluster. 

that the percolating cluster corresponds to the 0-£eld 


gested already in 15, 22 


It was sug- 
condensate 


4 

5 


The case I? = 4 is special because the gauge coupling is dimensionless. 

We do not discuss here subtleness which can arise from consideration of the coinciding points, xi 


in Eq. (^. 


X2, 
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in the classical approximation However, it is only the phenomenon of the hne tuning that 
makes the properties of the vacuum state non-trivial. Indeed, we have now two coexisting 
scales, a and Moreover, we will see that the properties of the percolating cluster 

differ radically from the properties of the hnite-size clusters. 


Lattice data 


We have already mentioned that the density of the monopoles in the percolating cluster 
scales, see (^. Recently, further measurements on the geometrical elements of the perco¬ 
lating cluster have been performed |T^. In more detail, the cluster consists of self-crossings 
and segments connecting the crossings. It was found that the average distance between the 
crossings measured along the trajectory scales, see Eq. (^). Violations of the scaling in 
Lfree negligible for all the lattices tested. 

One can also measure the average of the shortest, or Euclidean distance, (d), between 
the two crossings connected by segments. In this case violations of the scaling are more 
signihcant and, roughly, the data can be approximated as: 


{d)\/a K. (0.65 -h [{a^/a) — 0.25]) , 


(50) 


where a is the lattice spacing, a is the string tension for the SU{2) gluodynamics, and 
0.15 < a^/a < 0.35. 

Thus, we are confronted with the problem of interpreting the data (|[), (|^), (|50D . The 
hrst step is to reduce the set of data to a simple picture: normalized to a free-particle case 
the measurements correspond to an inhnite mass JS, ml ~ a“^. Which is at hrst sight 


a shocking observation defying everything that we said so far. Postponing discussion of the 
physics till the next subsection let us reiterate the reasons for such an interpretation of the 
data. 

For a free particle, one can measure its mass by comparing the Euclidean distance between 
two points with the length along the corresponding trajectory connecting the same points. 
Indeed, one obtains the distance L by differentiating the partition function (^) with respect 
to the chemical potential /i dehned in Eq. ([39|). On the other hand, the dependence of a free 
particle propagator on the chemical potential enters through the factor exp(—<^) where 
the mprop is as in Eq. ([T^). Differentiating with respect to fi the propagator D{mprop, d) of 
a free particle of mass mprop we get therefore: 


/r \ __ rl\ ~ *^/?-ee particle 

\-^free particle/ — Ail i_yi//tpT’op, tx i ~ / \ \ ’ 

a/i 2{mprop{a) ■ a) 


(51) 


where we neglected higher corrections in a and mprop (generally speaking, mprop also depends 
on a). 

Two particular cases are worth mentioning. First, if is in the physical units. 


"^prop ~ ^QCD then we would have 


(L) ~ 

a 


Literally, a model with a tachyonic mass was introduced first in ||22| and such a model would result in 
... ^-3 


Pmon 
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which is in blatant disagreement with the data (|50[) . Moreover if we assume cancellation of 
only the leading 1/a term in (|^ we would expect: 


which still cannot be reconciled with the data either. 

Finally, if ~ a“^ then 

(L) ~ {d) 

which does agree with the data (^) for a available as far as we neglect the correction linear 
in a. 

To summarize, we are coming to a paradox. Indeed, even if we accepted ~ 

to satisfy (L) ~ (d) this would not help since we would not be able then to explain the 
persistence of the physical scale in (|22|), (|50|) . 


B. Analogy to the Mofibauer effect 

Thus, we see that the hnite-size clusters and the percolating cluster exhibit different 
patterns of the monopole kinematics. In the former case we could neglect the monopole 
mass while in the latter case it is not possible at all. Looking for an analogy, we naturally 
come to the MoBbauer effect. Indeed, one could say that the effect is the difference in 
kinematics inherent to the decays of a free atom and of an atom belonging to a lattice 
in its ground state. Formally, the kinematics of the decay of the atom belonging to the 
lattice looks as if the atom had an inhnite mass. Thus, the analogy we are going to pursue 
is between the percolating monopole cluster and atoms in the ground state of a lattice of 
atoms. The reason for the analogy is that in the both cases we are dealing with ground 
states of systems with a gap. 

Let us hrst remind the reader the basic features of the MoBbauer phenomenon (for a 
review see, e.g., Ref. [^). One considers decay of an atom belonging to a lattice of atoms. 
Then there exists a probability Pq that the lattice is not excited at all, i.e. no phonons are 
emitted. Then the recoil momentum is transferred to the whole of the lattice and, clearly, 
this corresponds to an inhnite effective mass of the decaying atom. 

In the language of the quantum mechanics, the probability of the MoBbauer transition is 
determined by the following matrix element: 

Po = K^|exp(ip^ ■ (52) 

where |i) is the ground state of the lattice, '^atom is the position of the decaying atom and p.^ 
its recoil momentum. Then, there are two crucial quantities which determine the estimate 
of Po, namely the “naive” recoil energy R and the energy of the gap. For non relativistic 
kinematics, 

R = py2Matom , 

and the gap energy, ujgap- A rough estimate for Pq is as follows: 

Po ~ exp{-R/ujgap) . (53) 

Thus, when the recoil energy is much larger than the gap, R 3> ojgap Po tends to zero and 
kinematics is the same as for a free-atom decay. In this limit we would come back to the 
parton-like picture. 
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The lattice data on the monopoles indicate that it is the opposite limit, R -C oOgap-, which 
can be relevant to the monopoles (if the analogy is correct at all). 


C. Quantum propagator as Brownian motion 


The MoBbauer effect concerns kinematics in the Minkowski space. On the other hand, 
we are considering the monopole motions in the Euclidean space. Moreover, there are no 
decays of the monopoles of course and, at hrst sight, there is no physical quantity analogous 
to p.y. However, we will argue in this subsection that the famous Brownian picture for 
the quantum propagator provides a key to identify an analog to |p^| in the monopole case. 

The classical Brownian motion (see, e.g.^ Ref. 0) is the motion of a particle of mass 


Mbt which results from exchange of momenta with particles of a medium which are not 
observed and chaotic. Then the Brownian motion is a random walk with the step 


6p ■ 6 t 
Mbr 


where 6t is the average time between collisions and 6p is the average momentum transfer 
during the collision. 

In case of the quantum propagator one should rather think in terms of the Brownian 
motion of a particle attached to a string with a non-vanishing tension. This is the origin 
of the inertia, or of the mass term (for details see Sect. 4.1). As for a non-vanishing Ap it 
is now quantum in origin and can be estimated from the uncertainty principle. Since the 
monopole trajectory is measured on the scale a, generically Ap ~ a~^ and is parametrically 
large if a —0. 

This simplest estimate should be corrected, however, for the fact that one can introduce 
mass and kinetic energy only after coarse-graining see Sect. 4.1. Therefore, 


Ap 



(54) 


where n was introduced in Eq. (^7|) . Note that n does not depend on a and for practical 
estimates one can use 


^ (2-3) . 


D. Limit a —> 0: academic and realistic 

Coming back to our analogy with the MoBbauer effect, we conclude that in the limit 
a —> 0 we expect that the existence of the gap is not important and we can think in terms 
of the parton model. Indeed, in the relativistic kinematics, R ~ Ap ~ a~^ while the mass of 
“free monopole” niprop ~ and can be neglected compared to Ap at small a. Moreover, 
the gap is provided by the glueball mass, rrigi and does not depend on a at all. 

This conclusion could be foreseen: at short distances, i.e. in the limit a —0 one can 
neglect the effects of the binding of the monopoles. What is more surprising is that for 
realistic values of a we are in fact far from the academic limit a ^ 0. This is due to interplay 
of various numerical factors. Consider, for example, a = (3 GeV)“^ which corresponds to 
smallest a available. Then: 

m2 ~ ^ ~ 3 GeV2, (Ap)^ ~ 2 GeV^ (55) 

CL 
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where we used n ~ 5 (see Eq. (|5^) ) and the chemical potential (|39|) can be estimated with 
the help of Eq. (^) as p ~ (1.6 Thus the ordering of mprop and Ap is still reversed 

at the presently available lattices compared to the expectations in the a —0 limit. 

The onset of the parton model is expected when copious excitation of glueballs is favored. 
Since the glueball is relatively heavy, rrigi ~ 1.5 GeV the realistic values of Ap are far from 
satisfying the condition 

{Apf > {mprop + rugif . (56) 

Thus, we come to the following conclusion: 

The recoil-free geometrical picture of the percolating cluster is perfectly consistent with 
the guantum mechanics as far as the resolution (which is of order ^/na) is not too high. 
At better resolution we expect to loose both point-like monopoles and the geometrical picture 
itself. 

Actually, in the data of Ref. IQ] one can already see some hints on violations of the 1/a 
behaviour of the monopole held-theoretical mass M{a). Thus, we are inclined to correlate 
these hrst indications to appearance of the monopole structure [§] and the scaling violations 
in dav, see Eq. (|5I1|). Our prediction is that these indications become clearer at smaller a. 

One can roughly estimate the values of Ocrit at which the parton picture begins to prevail 
over the recoil-less kinematics as: 


^crit 


(10 ^ 30)GeV . 


(57) 


At such a one can also expect that the monopoles do not look point-like. Indeed, in the 
held theoretical language it is quite trivial that the dominance of the inelastic processes and 
appearance of the form-factor for an “elastic” process are determined by the same dynamics. 

Appearance of such a large mass scale as (|57|) might come as a surprise. It is worth 
noticing therefore that there exists accumulating, although indirect evidence on the relevance 
of large mass scales to the QGD phenomenology (see, e.g., Ref. and references therein). 
Now we see how such a scale can be build up on Aqcd, proceeding through the glueball 
mass and other numerical factors. 


VI. MONOPOLES AND GLUEBALLS 


Thus, the next crucial question seems to be the coupling of the monopoles to glueballs. 
In the most general way the monopole-glueball interaction can be studied by measuring the 
correlation function: 


K(x-y) = (0|p(i),p(i/)|0) , 


(58) 


where 

p{t) = 0+(x,t)(/(x,t) . 


Asymptotic at large distances is sensitive to a massive excitation [24 


lim (0|p(t),p(0)|0) = Cl +C 2 e (59) 


which is likely to be a glueball. 

As far as we approximate monopoles as point-like and conhne ourselves to an effective 
theory without gluons then correlators like (|5^) are given by a sum over closed monopole 
loops: 

Z{M) 


K{x - y) 


d3 id32 


31=32=0 


( 60 ) 
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where 

M{z) = ji8^{x - z)+ j28^{y - z) , 
and the partition fnnction is given by: 

” 4 ^/ <'*^0 / M(x{T))dT)U, ( 61 ) 

where {--^xo means averaging over all the closed paths. Farther details and dehnitions can 
be fonnd, e.g.^ in the review [^. What is important for us now is that the sum over the 
paths can be calculated numerically using the lattice data since the monopole trajectories 
are directly observed. 

The constant in Eq. is related to the Pmon'- 

C-l i^Pmon ■ ) ( 62 ) 

and is entirely determined by the percolating cluster. Determination of the glueball mass 
would require, on the other hand, sensitivity to the quantum corrections, or to the hnite-size 
clusters’^. 

It could quite well be so that keeping the closed paths alone would not be adequate to 
determine the glueball mass through (|5^) . This would be a signal that gluonic intermediate 
states are important and the theory is not unitary without inclusion of such states. Since 
the gluons are not detected directly on the lattice the pure gluonic intermediate states would 
look as breaks in the monopole loops. 

Thus, measurements of the correlator (^) would be very important to further understand 
the structure of the effective theory of the monopole interactions. 


VII. CONCLUSIONS 

We started with the picture according to which non-Abelian gluodynamics, when pro¬ 
jected onto the scalar-held theory via monopoles, corresponds to a hne tuned theory 0. The 
monopoles which we considered are dehned (“detected”) through the Maximal Abelian pro¬ 
jection. However, the mass scales which exhibit mass hierarchy are gauge independent. The 
scales are provided by the SU{2) invariant action per unit length of the monopole trajectory, 
on one hand, and by inverse “free path length” (see (^^), on the other. The observed hne 
tuning suggests the use of the geometrical language to describe the monopoles. 

In this paper we confronted the polymer approach to the ehective theories of the mag¬ 
netically charged helds with the lattice data on the monopole clusters. 

There are two types of clusters, namely, hnite-size and percolating. In the former case we 
have demonstrated that the length spectrum and the dependence of the size of the clusters 
on their length are well understood within the percolation theory. An alternative language 
is the simplest vacuum loop in the polymer language. 


7 


The percolating cluster alone at distances d ^ dav is representing a cluster of fractal dimension Dfr = 4, 
(see our discussion of the uncorrelated percolation in Sect. i) and cannot be sensitive to the glueball 


mass. 
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There is a striking difference, however, between the finite-size clusters and the percolating 
cluster. In the former case one thinks in terms of the vanishing mass of the monopoles (on 
the scale of the lattice spacing a): 

finite—size ~ 0 ■ Cl . (63) 


This picture is conhrmed by evaluation of the critical exponents just mentioned. On the 
other hand, the straightforward interpretation of the geometry of the percolating cluster 
leads to 
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percolating ~ ® 


-2 


(64) 


We have argued that the apparent discrepancy between the effective monopole masses 
extracted in the two ways is naturally resolved if one takes into account that the hnite-size 
clusters correspond to quantum corrections while the percolating cluster corresponds to the 
ground state. The quantum corrections reveal then presence of a massless (on the scale of 
a) excitation. In case of the ground state, we observe an analog of the MoBbauer effect. The 
onset of a parton-like description is delayed numerically by some factors, most notably by a 
relatively high glueball mass rrigi where rrigi plays the role of the gap. 

We also predict that at much smaller a (see (|^) one should see the convergence of (|B1D 
to (^. 

Thus, hrst applications of the polymer picture to the monopole clusters turn successful. 
However, justihcation for such an approach remains at present pure phenomenological in 
nature and relies on the observation of the hne tuning. Any further theoretical development 
is therefore to be checked by measurements on the lattice. In particular, studies at smaller 
a, both theoretical and numerical, are desirable. 
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